Partition Functions of Strongly Correlated Electron Systems as "Fermionants" 
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We introduce a new mathematical object, the "fermionant" Fermjv(G), of type N of an n x n matrix G. It 
represents certain n-point functions involving TV species of free fermions. When N = 1, the fermionant reduces 
to the determinant. The partition function of the repulsive Hubbard model, of geometrically frustrated quantum 
antiferromagnets, and of Kondo lattice models can be expressed as fermionants of type N — 2, which naturally 
incorporates infinite on-site repulsion. A computation of the fermionant in polynomial time would solve many 
interesting fermion sign problems. 
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The Pauli exclusion principle has profound physical con- 
sequences, ranging from the stability of condensed matter 
to the structure of atoms, nuclei, as well as of protons and 
neutrons. In quantum mechanics, the Pauli principle man- 
ifests itself through the requirement that the wave function 
must be antisymmetric under the exchange of two fermions. 
In many-body quantum mechanics this is incorporated by 
the fact that fermion creation and annihilation operators anti- 
commute with one another. When a fermionic path integral 
is set up in an occupation number basis, the resulting world- 
lines of identical fermions form closed loops wrapping around 
the periodic Euclidean time direction. In relativistic theories, 
these worldlines can also form closed loops which incorporate 
the physics of pair creation and annihilation. The number of 
closed loops, which corresponds to the number of cycles in a 
permutation of fermion positions, determines the fermion per- 
mutation sign of a configuration of worldlines. Alternatively, 
fermion path integrals can be constructed in a basis of coher- 
ent states using anti-commuting Grassmann variables. Grass- 
mann variables can be viewed as a different way of bookkeep- 
ing of the fermion worldlines. Indeed, Feynman diagrams are 
well-known representations of these worldlines and the nega- 
tive sign associated with each fermion loop is nothing but the 
fermion permutation sign. 

While worldlines of bosonic particles can be given a clas- 
sical interpretation, the same is not true for fermions due to 
the fermion permutation sign, which is an important quantum 
mechanical phase that often makes the weight of a configura- 
tion of fermion worldlines negative. This makes it impossi- 
ble to interpret this weight as a probability in a Monte Carlo 
simulation based on importance sampling. Fortunately, it is 
possible to sum over all fermion worldlines when fermions 
are not interacting with one another. For this reason, the 
traditional aim of Monte Carlo approaches to fermions is to 
convert an interacting into a non-interacting problem, by in- 
troducing auxiliary fields which incorporate the interactions. 
Once the fermions are freely moving in the auxiliary field, 
one can sum over all fermion worldlines which leaves be- 
hind a fermion determinant as the weight of each auxiliary 
field configuration. In cases where this determinant is posi- 
tive, the sign problem is considered solved and a Monte Carlo 



algorithm can be devised to sample the auxiliary field. This 
has been the method of choice for almost all strongly corre- 
lated fermion problems until recently. However, in general 
there is no guarantee that the fermion determinant is posi- 
tive. The physics of the auxiliary field plays an important role 
in determining its sign. There are many interesting quantum 
systems, ranging from high-temperature superconductors, ge- 
ometrically frustrated antiferromagnets, and Kondo lattices, 
to systems of quarks and gluons at non-zero baryon density, 
where indeed the fermion determinant can be negative or com- 
plex. We refer to these as difficult sign problems that remain 
unsolved. If only the magnitude of the weight of a configura- 
tion is used for importance sampling while its sign or complex 
phase is included in measured observables, the error to signal 
ratio and thus the required computational effort increases ex- 
ponentially with the space-time volume. 

Over the last decade new ways to sum over fermion world- 
line configurations have emerged which allow us to write the 
partition function as a sum over novel types of configurations 
with non-negative weights. In the meron-cluster approach, 
fermion worldlines are divided into independent clusters. A 
summation over the fermionic degrees of freedom within a 
cluster produces a non-negative weight for each cluster con- 
figuration lll|2(]. More recently, this idea was extended to the 
notion of fermion bags, where again a summation of fermion 
worldlines within a bag yields a non-negative weight for each 
fermion bag [3]. In yet another development, in certain mod- 
els a sum over all Feynman diagrams of a given order in the 
perturbative expansion can be performed and yields a non- 
negative weight for these diagrams. This leads to the diagram- 
matic determinantal Monte Carlo method ||4j,|5[]. While these 
new approaches allow us to solve some sign problems, which 
seemed unsolvable with the traditional auxiliary field method, 
none of these new developments have yet solved some of the 
difficult sign problems mentioned above. 

In this letter we further investigate the nature of these diffi- 
cult fermion sign problems by showing that the partition func- 
tion of some strongly correlated electron systems, including 
the infinite U Hubbard model away from half filling, the quan- 
tum Heisenberg model for geometrically frustrated antiferro- 
magnets, as well as the Kondo lattice system, can be expressed 
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in terms of a new mathematical object — which we name the 
"fermionant" — which is a generalization of the determinant 
of a matrix. The exact evaluation of the partition function is 
a sum over exponentially many states in the Hilbert space and 
so is generically impossible in polynomial time. Remarkably, 
for free fermions, since the partition function is a determinant 
of an n x n matrix it can be computed with an effort propor- 
tional to n 3 . Our results imply that some interesting strongly 
correlated electron systems share certain features with free 
fermions and that in these systems the determinant is replaced 
by a fermionant. However, it is not clear whether the fermio- 
nant can be computed in polynomial time. Thus, we think 
that perhaps the computation of a fermionant is at the heart 
of the difficult sign problems. Fermionants naturally arise in 
systems in which different fermion species (such as spin up 
and spin down) are not allowed to simultaneously occupy the 
same lattice point. Interestingly, fermionants also describes a 
class of n-point functions in theories of free fermions. 

Let G be an arbitrary real- or complex-valued n x n matrix. 
The fermionant of type N of the matrix G is defined as 

Ferm^G) = (-1)" £(-JV)"Mg 1(t(1) G 2(t(2) . . . G na(n) , 

(1) 

where the sum extends over all permutations a of the indices 
1,2,. ..,n, u(a) is the number of cycles in the permutation, 
and cr(i) is the permutation partner of the index i. Based on 
this definition, it is clear that the fermionant of type N = 1 is 
simply the determinant. 

Let us consider N species of fermions moving on a space- 
time lattice with sites labeled by the index x = 1, 2, . . . , V, 
where V is the space-time volume. We use Vx,<* and i() x a 
to denote the Grassmann variables associated with the species 
a = 1, 2, . . . , N of the fermion at the site x. We assume that 
the fermions hop on the space-time lattice as non-interacting 
particles, based on the [/(iV)-invariant Euclidean action 

N 

S$, <p] = J2 D«,« M ^,«. ( 2 ) 

q = 1 x,y 

which results in the partition function 

Z= f #*,«] eM-S®,i>}). (3) 

J x,a 

We now define the n-point function 

C(x 1 ,X 2 ,...,X n ) = [ J|[# x , a #x,a] 

J x,a 

x^ Xi ip Xl ... ^ x J) Xn exp(-5[^,^]). (4) 

In this expression ip x .ip Xi = ^2 a ip x . aV'xi.a contains a sum- 
mation over the species index which has been suppressed 
for convenience. Using Feynman rules, one can show that 
C(x\,X2, ■ ■ ■ ,x n ) = Fermjv(G), where G is the n x n ma- 
trix composed of elements connecting sites Xi and Xj given by 



G XiXj = (M _1 ) XiXj , where M _1 is the free fermion propa- 
gator. 

It is well-known that free fermion partition functions can be 
written as determinants of V x V matrices. It is also known 
that partition functions of interacting theories are sums over 
different n-point functions of a free theory. This can easily be 
seen in perturbation theory by expanding in powers of the cou- 
pling and then evaluating each term as a correlation function 
in a free theory. For some familiar strongly correlated elec- 
tron models, we now show that the partition function can be 
expressed as a single U-point function of a free theory, which 
turns out to be a fermionant of type iV = 2ofaV r xV r matrix. 

Before discussing the partition functions, it may be use- 
ful to recall the coherent state path integral representation of 
fermions using Grassmann variables. Let ct a and cs, a rep- 
resent the creation and annihilation operators for electrons of 
spin a =t, I at the spatial lattice site x. Introducing Grass- 
mann variables ip StCt and ip S a , an electron coherent state at 
the site x is defined as 

|Vx) = exp(V^t4t + ^4i)l°)> 

(4>s\ = (0|exp(c st VxT + c £4.V'xu)- ( 5 ) 

It is easy to verify that ($ s \ip s } = exp(^ )f ^ )t + ^ Sii ^s,i.) 
and that the completeness relation is given by 

|0)(0| ^| 0)(0 | c . t + 4j )(0|c^ + 4 t 4lO)(0|c^ C;?t 

J a 

where we have defined ipgips = V'x.t^.t + V'xuV'x.j.- If one 
wants to remove the doubly occupied state from the Hilbert 
space, which implies an infinitely strong repulsion between 
the electrons of spin up and spin down, this is easily accom- 
plished using 

|0)(0|+4 t |0)(0|c st + 4j0)(0|c^ = 

If one also wants to forbid the completely empty state, one 
can accomplish this using the relation 

ci t \0}(0\c Sf + 4j0)(0|c^ = 

n[di>x, a di>x,a] (~1Px^s) IV'xXV'xI-C 8 ) 
J a 

Thus, by introducing different types of completeness rela- 
tions, one can reduce the Hilbert space of the non-interacting 
theory and thereby introduce interactions. The partition func- 
tions of the resulting constrained theories are then equivalent 
to a single correlation function of a non-interacting theory. 
Based on this insight we now consider specific partition func- 
tions. 
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Let us consider the lattice action of free non-relativistic 
electrons hopping on a space-time lattice given by 

^Hubbard = — ) ^ ^ x ,r +\ ^ x ,t ,a ~ 'x ,t ;ai- 'x \r ,a 

x,T,a \ 

+ Ste^ ^ (l>*,T+l, a *l>g+l,r, a + ^x+lr+l,a^r, a ) > (9) 
i ' 

where (x, r) is a Euclidean space-time lattice site and ipx,r,a 
and ip s T a each represent Grassmann variables associated 
with the electron of spin a at the corresponding space-time 
point. The parameter e represents the lattice spacing in Eu- 
clidean time, while t and /i are the hopping parameter and the 
chemical potential, respectively. To obtain the Hamiltonian 
version of the problem, one needs to take the e — > limit, 
keeping eL t = 1 /T fixed, where T is the temperature and L t 
is the number of lattice points in Euclidean time. Based on the 
previous discussion of the completeness relations of coherent 
states, it is easy to show that the V^-point function 

Z = J [dip dip] II ^ ~ V'a.T^.T^ exp(-S'Hubbard) (10) 

is the partition function of the Hubbard model with infinite re- 
pulsion U between spin up and spin down electrons. Note that 
we have used eq. (Q just to project out the doubly occupied 
state from a free electron theory. 

At half-filling, the infinite U Hubbard model naturally re- 
duces to a spin | quantum antiferromagnet, provided contri- 
butions proportional to i 2 are maintained. In the coherent state 
path integral this is accomplished if one uses the action 



'magnet 



^Y,(^S,r+l, a ^3+lr, a +^x+lr + l^,r, a ) >, (H) 
i ' 

instead of Sjjubbard in eq. ( [Tol l. Note that the spatial lattice 
remains unspecified and could hence be either bi-partite or 
geometrically frustrated. 

Finally, a model relevant for the physics of heavy fermions 
is the so-called Kondo lattice model [6]. In this model a bath 
of free electrons interacts with a lattice of localized spins or 
magnetic moments. We construct this model by adding two 
free fermion actions. The free electron action is given by 



'electron 



E 



i ' 

where double occupancy is now allowed. Additional local- 
ized spins are described by other Grassmann fields Xr T a an d 



Xfi T a associated with the lattice sites R on which impurities 
are located. We denote the space-time volume occupied by 
the impurity sites by Vr. The free impurity action is given by 



impurity 



R,r,a 

VeJ (%A,T+l,J>H,r, a + ^a,T+l,aXa, r> a)|- ^ 

The partition function of the Kondo lattice model can now be 
written as the Vr -point function 

z = J [##] [dxd x ] n ( -Xr, t Xr iT ) 



cxp(-5 c lectron ^impurity) 



(14) 



Since the localized impurities can either be spin up or spin 
down, we now eliminate the doubly empty and doubly occu- 
pied states only of the impurity spins. 

In the three cases discussed above, the partition function is a 
V^-point (or in the case of the Kondo lattice model a Vjj-point) 
function of a free theory, which is given by a fermionant of 
type N = 2, 



Z = Z Ferm 2 (G), 



(15) 



of a V x V (or a Vr x Vr) matrix G — A/~\ which 
represents the fermion propagator of the corresponding free 
fermion lattice action written in the form of eq. (O. The factor 
Z = Det(M) 2 is the partition function of the free lattice ac- 
tion. Some observables, such as certain correlation functions 
of fermion bilinears, can also be expressed as fermionants. 

By rewriting a well-known NP-complete problem as a spin 
model with a sign problem, it has been shown that some sign 
problems are as hard as NP-complete problems 0], These 
problems belong to the hardest problems in the complexity 
class NP [8]. Since it is generally expected that NP ^ P, this 
result implies that a generally applicable Monte Carlo method 
that solves all sign problems with autocorrelation times that 
scale at most polynomially with the system size should not 
exist. This does not exclude that specific sign problems can 
still be solved. It is important to understand whether the phys- 
ically relevant difficult sign problems mentioned before are 
also NP-complete, or whether one may hope for a solution in 
polynomial time. If we consider a sign problem solved when 
the partition function can be expressed as a sum over config- 
urations with non-negative weights, such that the weight of 
each configuration can be computed in a time that grows at 
most polynomially with the system size, then our work shows 
that the solution of difficult sign problems is closely related to 
the computational complexity of fermionants. 

Littlewood and Richardson have introduced the concept of 
the immanant 



Immr(G) = V Xt{v)Gi<j(i)G 2 o{2) ■ ■ ■ G T 



r(n)> 



(16) 
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as a generalization of the determinant of an n x n matrix G |@1 ■ 
Here xr(&) is the character of the permutation group element 
(j in the representation T (characterized by a Young tableau 
with n boxes). When one considers the totally antisymmetric 
representation (associated with the Young tableau consisting 
of a single column of n boxes), the immanant reduces to the 
determinant. Similarly, when one considers the totally sym- 
metric representation (associated with the Young tableau con- 
sisting of a single row of n boxes), the immanant reduces to 
the permanent Perm r (G) = J2 a Gia(i)G 2 a(2) ■ ■ ■ Gna(n)- 
The evaluation of the permanent is even more difficult than 
solving NP-complete problems ifioll . and the computational 
complexity of a general immanent is similar to the one of the 
permanent [11]. Remarkably, thanks to its invariance proper- 
ties, the closely related determinant, whose naive computation 
would require an effort proportional to nl, can be calculated 
in just n 3, steps. 

Just like the permanent and other immanants, the fermio- 
nant of a matrix is a generalization of the determinant. The 
fermionant weighs the contribution of a permutation a with 
the factor (— N)"^, where v{a) is the number of cycles of 
a. Just like the character xr(c), the number of cycles v(a) is 
the same for all members of a conjugacy class. The characters 
form a complete orthogonal system. Consequently, one can 
expand 



(-AT (CT) =E c rWxr(a) 



(17) 



such that the fermionant can be expressed as a linear combi- 
nation of immanants 

Ferm N (G) = (-l) n ^ c r (N)Imm r (G). (18) 



In these systems electrons of spin up and spin down are not 
allowed to simultaneously occupy certain lattice points. Re- 
markably, the partition function of such systems is naturally 
described by a new mathematical object — the fermionant of 
a large matrix — which is a simple generalization of the de- 
terminant. Interestingly, the fermionant also represents a par- 
ticular rt-point function of a system of free fermions. In par- 
ticular, fermionants naturally appear when one sums over all 
Feynman diagrams of a given order in the perturbative expan- 
sion of the repulsive Hubbard model. The fermionant is math- 
ematically related to the immanents of a matrix, for whose 
computation — in contrast to the determinant — no polyno- 
mial time algorithm is known. If a polynomial time algorithm 
for computing the fermionant would exist, several physically 
relevant fermion sign problems could be solved. While this 
may not be the case, the relation between the fermionant and 
the partition function of these strongly correlated electron sys- 
tems sheds new light on the computational complexity of the 
corresponding fermion sign problems. 
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While this relation is not useful for numerical purposes, it 
shows that the fermionant is mathematically intimately related 
to the immanants. Although no polynomial time algorithm is 
known for the evaluation of a general immanant, this does not 
necessarily mean that the search for a polynomial time algo- 
rithm to evaluate the fermionant is doomed from the outset. 
In particular, a careful examination of its invariance proper- 
ties seems worthwhile. Even if the evaluation of a general 
fermionant may require a non-polynomial computational ef- 
fort, some fermionants can still be sampled stochastically in 
polynomial time. For example, non-frustrated quantum an- 
tiferromagnets, whose partition functions indeed are fermio- 
nants, can be simulated with the very efficient loop-cluster al- 
gorithm B12lll3ll . Since it arises in several important physical 
systems, a careful study of the mathematical properties of the 
fermionant is highly desirable. 

In conclusion, we have investigated systems of strongly in- 
teracting fermions including the infinite U Hubbard model, 
quantum antiferromagnets on arbitrary bi-partite or geomet- 
rically frustrated lattices, as well as Kondo lattice models. 
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